Non-equilibrium 2-Channel Kondo model for quantum dots 
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Wc find that, under certain condition, a quantum dot with odd number of electron and coupled 
to two leads can be described by a non-equilibrium 2-channel Kondo model, when the two leads 
has a large voltage bias between them. The model is exactly soluble and can be mapped into a free 
fcrmion system, even in the presence of external magnetic field and other relevant perturbations. All 
(dynamical) correlation functions can be calculated. The fixed point of the 2-channel Kondo model 
is a non-equilibrium fixed point and is different from the usual 1-channel and 2-channel Kondo fixed 
point. 

PACS numbers: 73.23.-b, 71.10.Pm 
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In a recent experiment, Goldhaber-Gordon etal jy and 
Cronenwett etal |2) observed a Kondo effect in quantum 
dot systems predicted by theories. |^-[ic|] The Kondo ef- 
fect in quantum dot systems is very interesting since 
many parameters can be adjusted, which allow us to 
probe many different regime/aspects of the Kondo ef- 
fects. 

In this paper we are going to study the Kondo effect 
in quantum dot systems in the limit where the charge 
fluctuations on the dot can be ignored, and when there 
is a large voltage bias across the two leads. Under cer- 
tain conditions as will be indicated bellow, the properties 
of the system is described by a 2-channel Kondo model, 
with non-equilibrium effect included as a perturbation. 
This 2-channel Kondo model can be solved exactly us- 
ing a Major ana fermion approach even in the presence 
of some relevant perturbations, such as finite tempera- 
tures, finite magnetic fields, non-equilibrium tunneling 
between leads, and finite changes in the relative strength 
of the coupling between the two leads and the dot. Many 
properties of the system (including dynamical proper- 
ties) can be calculated exactly. The fixed point of the 2- 
channel Kondo model studied here is a non-equilibrium 
fixed point which is different from the usual 2-channcl 
Kondo fixed point. For example the dot-spin correlation 
behaves as t~ 2 here instead of t~ Y for the usual 2-channel 
Kondo fixed point. 

We start with the following model Lagrangian for the 
quantum dot 



C 



iiph( d t + vd x + i-eVrf^ipja - h ■ SS(x) 
-2XS(x)S ■ (u*4 a + v*4 a )?f-(ml> w + whp) (1) 



where V is the voltage difference between the first and 
the second lead, and h is the coupling of external mag- 
netic field to the dot-spin. Note that the spin-up and the 
spin-down Fermi surfaces in the leads remains to have 
the same energy even in the presence of external mag- 
netic field, thus there is no coupling between the external 
magnetic field and the lead-spins in our model. Am and 



Xv describe the strength of coupling between the dot and 
the two leads. 

After redefine the ip field in Eq. ([[]): tpi a — > 
ipj a , the above Lagrangian can be rewritten as 



£ = £ 



Knd 



(2) 



C-Knd = WiaiPt + vd x )lp ja 



-2X5(x)S ■ (MVJa^i/s + MVL^ 2/3 ) (3) 

Ct - -2X t 6(x)S ■ (e iey Vw4 a ^V 2 /3 + h.c.) (4) 
-h ■ S6{x) (5) 



C H 



where A t = A. We see that the finite bias is described by 
a time dependent term after the mapping. At is the am- 
plitude of the spin-flip tunneling between the two leads, 
and A is the spin-exchange coupling between the dot-spin 
and the lead-spins. 
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FIG. 1. The tunneling can be reduced by a change in the 
density of state. 

The above system is a 1-channel Kondo system when 
eV = since A f = A. However, when eV is large, the 
spin-exchange term and the spin-flip tunneling term be- 
come very different. If one can make At <^ A in the large 
eV limit, then the system becomes a 2-channel Kondo 
problem. One way to reach this limit is to use the energy 
dependence of the density of states. For example, one 
may use n-type semiconductor as one lead and p-type as 
the other lead {ie put the quantum dot in the depletion 
layer of a diode). When eV > Ep, we have A(<A since 
the tunneling is blocked (see Fig. |IJ). Certainly, to see 
the Kondo effect, we need Ep to be less then the level 
spacing in the dot. 
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Even when A = At at high energies, we can effectively 
have A > A t at low energies. This is because at energy 
scales above eV, the system is a 1-channel Kondo sys- 
tem, and A and At flow together to larger values as we 
lower the energy scale. When energy scales are less then 
eV, At can no longer flow. If the value of At is small 
enough at eV (ie if eV 3> Tk where Tk is the 1-channcl 
Kondo temperature for eV — 0), the system behave like 
a 2-channel Kondo system and A can continue to flow to 
even larger values as we lower the energy scales further. 
Thus at low energies, we effectively have A t > A. How- 
ever, since A and At flow as ln(-E'), it may be difficult to 
achieve A 3> A t at low energies. 

When A t <C A (at low energies), we may first set 
At = 0, and study the 2-channel Kondo problem de- 
scribed by Cxnd- The 2-channel Kondo system CkuJ, 
has a 2-channel Kondo fixed point, which can be reached 
by tuning the relative strengths of the couplings between 



the dot and the two leads to 



, and lowering the 



temperature below a 2-channcl Kondo temperature T 2 k- 
Then, we can include back the £t term to study its effect 
on the 2-channel Kondo fixed point. As we will see later 
that, at the fixed point, the effects of the Ct term (as well 
as many other terms, such as the coupling of dot-spin to 
the external magnetic field) can be calculated exactly. 

In the following we are going to introduce a Majo- 
rana fermion approach (which is a combination of the 
current algebra approach in Ref. [lT|l2]| a nd the Majo- 
rana fermions approaches in Ref. pl[ l3| , |l5| ]) to solve the 
2-channel Kondo fixed point. The Majorana fermion ap- 
proach used here has a full spin rotation symmetry, and 
is similar to the one used in Ref. fl6| ]. The 2-channel 
Kondo Hamiltonian is given by 

H = H f + H K = iv J dx (i>\ a d x ^ ia + AJ s (0) • s) (6) 

where Hf is the free fermion Hamiltonian. The above 
2-channel Kondo system flows to a fixed point with cou- 
pling A = Ao = irv. In the rest of the paper we are going 
to only consider physics at energy scale much lower then 
the Kondo temperature T 2 k and concentrate on the fixed 
point Hamiltonian with A = Ao. 

Introducing the charge U c (l), spin SU S (2), and flavor 
S'[/ C h(2) densities (we will call the SU c h(2) quantum num- 
bers as flavors): 



t ' T "/ 3 ,/,._ T A _ „,,t T ij 



we can rewrite the 2-channel Kondo Hamiltonian as 



-J 



H K = XS(x)S ■ J s 



(8) 



Introducing J s (x) = J s {x) + S(x)S which satisfies the 
same algebra as J s , |ll[] we can rewrite the above Hamil- 
tonian as 



#2CK - —J 2 + —J s + —{Jf? + const. (9) 



irv ~ 2 

Y 



irv 

T 



if A = Ao. In terms of J s , the 2-channel Kondo Hamilto- 
nian at the fixed point A = Ao takes the same form as a 
free fermion Hamiltonian. 

In the following, we are going to use Majorana fermions 
to describe our 2-channel Kondo system. Introduce eight 
Majorana fermions \^2Keip c , \/2Im?p c , r/ A |a=i,2,3 and 
Vs\a=i,2,3, |l3| where ip c carries y/2 unit of the U c (l) 
charges, ?7^|a=i,2,3 form a flavor triplet and ?7°| Q =i,2,3 
a spin triplet. The free Hamiltonian Hf can be rewritten 
as 



H f =ivjdx (v40*V>c + lvfd x Vf + \ri a s d x ri1 



(10) 



The free theory is described by the SU S (2) 2 x S , [/ ch (2) 2 x 
U c (l) KM algebra. In terms of the Majorana fermions, 
the currents J c , J s , and J A take the form J c — v / 2'0c'0c, 

ja = e ab Cjl a v bJ 2: an d jj = e ABC tfrf /2. 

Let us concentrate on the spin part described by 



H s 



IV 



(11) 



Assume the system is finite —L/2 < x < L/2, and 77" 
satisfy the anti-periodic boundary condition: r]^(x) — 
—rj^{x + L). In this case the ground state |0) is a spin 
singlet. The total Hilbert space contains two sectors: 
Tt s = Hq © Hf . Hq containing even number of fermions 
is the space generated from |0) by applying any number 
of the J s operators, while H{ containing odd number of 
fermions is generated by applying any number of the J s 
operators and one operators. 

Now let us add the Kondo term Hk at the critical cou- 
pling Ao- The spin part is still described by the S , C/ S (2) 2 
KM algebra in terms of the new SU S (2) current J s . No- 
tice that Hk conserves the even-odd fermion numbers. 
Thus the total Hilbert space still contains two sectors: 
H s = r H.\i 2 © ^i/2> where (or W*y 2 ) carries even- 

number (or odd- number) of fermions. The system can 
be solved using KM algebra. The ground states in each 
sector carry spin 1/2 due to the added dot-spin. All the 
states in each sector are generated by applying J s 's on 
the corresponding ground states and the states in Ti.^ 2 
and H*/2 have a one-to-one correspondence. We see that 
the Kondo coupling induces the following simple changes 
in the Hilbert spaces: p2] | 



n -> n 1 



/2> 



HI 



n_/Sf 

H-1/2 



(12) 



We can introduce new Majorana fermions to represent 
the new SU S (2) current: J° = ^e abc fi^fi h s . f)" carry 
odd fermion numbers and map between even- and odd- 
fermion sector. When x ^ 0, 77° (x) = sgn(x)r]® (x) . 
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Note that the total Hilbert space at the fixed point has 
a structure H s — TCy 2 ® {|0), |1)}. We can introduce a 
fermion operator eft that map the even-fermion state |0) 
to the odd- fermion state |1): |1) = d'\0). Note that d 
and d) commute with J s , and do not appear in the fixed 
point Hamiltonian (which contains only J s ). Thus d is 
a free fermion operator which carries no energy. Using 
the commutation relation (which actually defines a spin 
triplet primary field) 

[J2(x),S b ] = ie abc 8 (x)S c (13) 

and the fact that S carry even-fermion numbers, we find 
that the dot-spin operator S is given by 

S = ^Wj.(0) (14) 

where b = d + d' . The normalization coefficient is ob- 
tained by noticing that {^(0), ^(0)} = 6 ab 6(0) and 
(5(0) = 1/a if we choose a finite short distance cut-off 
a. With the fixed point Hamiltonian 

H s + H K = iv J dx \f, a s d x r, a s (15) 

we can now easily calculate the S correlation. 

Now let us add back the non-equilibrium tunneling 
term C t and include the channel asymmetry term 

8-^1(0)4*^(0) (16) 

Both terms are spin singlet and flavor triplet. Thus the 
leading contribution to those operators at the 2-channel 
Kondo fixed point comes form the flavor triplet primary 
fields which have dimension 1/2. To obtain the Majorana 
fermion representation of the flavor triplet primary fields, 
let us analyze the Hilbert space of the spin and the flavor 
sectors (with zero-charge). Without the Kondo coupling, 
the Hilbert space has a form H sf = H s ® n{ © H{ ® n{ , 
where H.q ± are defined above for the spin sector, and 
Til 1 are defined similarly for the flavor sector. The sec- 
tor H s Q (g)H f (or H{®H{) is generated by J s a and Jf from 

|0) (or rfirif\a) ~ ^Vy^sl ))- After we include the 
Kondo coupling, according to Eq. ([l2]) , the Hilbert space 
becomes H sf = H s 1/2 ®n{® H{', 2 <g> H{. The above 
structure of the Hilbert space tells us that the spin triplet 
Majorana fermion operator jj" are unphysical since they 
map a state in "H\i 2 ® T~(-o to a state in Hf, 2 Hq which 
is outside of our physical Hilbert space. Similarly the 
flavor triplet Majorana fermion operator are also un- 
physical. The only physical spin triplet primary fields 
are S and the only physical flavor triplet primary fields 
are br)f- We note that the flavor triplet primary fields 
are quadratic in the fermion operators. Our model for 
the fixed point remains to be a free fermion theory even 
after we include the relevant perturbations represented 



by brjf- This allows us to study exactly the effect of Ct 
term and the channel asymmetry term on the 2-channel 
Kondo fixed point. 

More generally, we find the following exactly soluble 



model 

H = H 2 ck + H a + H' a 

+H t +H' t + H h + H' h + H / fs (17) 

#2CK = -iv J dx + VfY^f + C yC) (18) 

H a =i la bn}(0) (19) 

K = Tai<? AB rif(P)rif(0) (20) 

H t = i^= lt e ieVt hp){0) + h.c. (21) 

H' t = z-^=e^V}(0) (7^(0) + 7 U 3 (0)) + h.c. (22) 

H h =ibrf a .fj.{0) (23) 

H' h = j' s a ie abc r) b s (0K(0) (24) 

H' fs = i^tvf(0K(0) (25) 

where V/ = (vj ~ i? ?/)/%/2 and 7 S = ^h. This 



Hamiltonian describes what we call the non-equilibrium 
2-channel Kondo model. H 2 qk in H describes the 2- 
channel Kondo fixed point. All other terms represent 
relevant or marginal perturbations around the 2-channel 
Kondo fixed point. H a and H' a are induced by the asym- 
metry in the coupling between the channel spins and the 
dot-spin. Ht and H' t are the non-equilibrium tunneling 
terms caused by the finite voltage bias between the two 
channels. H t corresponds to the spin- flip tunneling de- 
scribed by Ct and H[ corresponds to the direct tunneling 
described by ^LV^a + h.c. and i/>l a a^^20 + h.c. H h is 
the coupling between a magnetic field h and the dot-spin 
S. 

To calculate the physical properties of the above sys- 
tem we also need to specify the "boundary" condition, 
ie how incoming fermions are distributed. The bound- 
ary condition here is that all income branches contain no 
Majorana fermions, ie all k > levels are empty. 

Next we simplify the Hamiltonian using the transfor- 
mation ipf(x,t) — > e leV (t+ v x ^ipj(x,t) which changes 
the eV in H to zero. However, the transformation 
also changes the boundary condition for the incoming 
branches: ipf is now filled up to energy eV, but the 
k > levels for all other Majorana fermions remain to 
be empty. 

For simplicity, let us drop the sub-leading terms H n s. 
Among the eight Majorana fermions in H (now with 
eV = 0) only one linear combination couples to b: 

f d x 
H b = dx (-iv)ri b —T] b + ijbr] b (0) (26) 

where 
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% = 7 1 [7 '?/+Re(7t)7 ? }+Im(7t)^+7 S ^ S ] (27) 
- = \ '7, - |7t| 2 +7 2 (28) 



V / 7a + l7tP+7s 

The equation of motion are |l7| 

dtr? 6 (x, t) = -vd x r) b (x, t) - jb(t)5(x) (29) 

d t b(t) = -2 7 %(0,f) = -7(%(0 + ,t) + %(0",0) (30) 

If we expand ry(x) = ^2 k rj_(k)e tkx for x < and rj(x) = 
Ylk^+i^) 6 for x > 0, then the scattering is simply an 
energy dependent phase shift: r/ + (k) = e l ^ k rj^(k) where 



iv k 



'j 2 k 



■T 



(31) 



This allows us to determine how the eight Majorana 
fermions, fa, ry 8 ) = (v^R.ei/'c, */2bm() c , rjf , 77" ), scat- 
ter across the d level at x = 0: 



7 1I+ {k) = M IJ T 1 j„{k), M=l + (e 



i<Pk 



1 



7/7J 



(32) 



where ( 7 i,...,7 8 ) = (0, 0, Re(7 t ), Imfa), 7o , -fl, 7^, 7 f), 
and Tjj±{k) are the Fourier components of 77/ (x) for x > 
and x < 0. The large phase shift <r3 fe = 7r at fc = 
indicates a resonance at the Fermi level. This setup al- 
lows us to easily calculate the total tunneling current JTt]] 

It =evf& (^J + (fc)^ /+ (fe)-^J_(*)^-(*)). We find 
that 

«l7t| 2 



2^(M 



-7a 2 



7 2)2 + u 4 fc 2 



(l7 t | 2 + 7 2 +7 a 2 )K^ + e^)-n(«/c- e y)) (33) [1] 



where 11(e) — l/(e c / T +l). Unfortunately I t having a very 
smooth dependence on V cannot reveal the resonance at 
the Fermi surface when eV is large. One need to use the 
noise in I t near w = eV to probe the resonance. E3] 

The effect of H' is to generate additional energy in- 
dependent rotations among r\i 's. It changes M to M = 
1 + (e*^ fc — 1)7/77, wnere 7 ± are t w0 un it vectors. 

The equation of motion for b also implies that b = 
-i J H^(l + e^ fc )?7_(fc) which satisfies b 2 = 1. This 
expression allows us to calculate the response function 
iD f (t) = Q(t)([F 3 (t),F 3 {0)}) where F 3 is the chan- 
nel asymmetry operator F 3 = ibn 3 (which comes from 



Eq. (|I6|)). We find that (assuming 7 a 
eV — > 00 for simplicity) 



Is 



ImDfM = 

n ' 2ttv 



dv 



9 _ 
7 v 



D f (u) = ~J dv 



v 2 + 7 4 u 
ImD / (^) 
j - v + i0+ 



7 (n(v - u) 



and 

■ n(w — u)) 
(34) 



Assume 7 a depends on a gate voltage V g which control 
the strength of coupling between the dot and the leads. 
Changing V g will generate a term AV g ! ^f L F 3 in Hamil- 

tonian. Thus Q = ^f-F 3 is the charge operator that 
couples to V g . The induced charge is 



(35) 



We see that C — {jjy^J Re£>/(w) is the capacitance, 

and a = u> (^§y^j ImDf(cj) is the conductance in paral- 
lel with the capacitance. We find that, for u> > 7 2 u _1 , 
C r~j ln(T 2 ^/tj) and a cx u>, while for u> < 7 2 u _1 , 
C ~ ln(wT2i<"/7 2 ) and er oc oj 2 . Measuring C and a by 
applying an AC component in the gate voltage V g may 
allow us to probe the resonance associated with the fixed 
point. The behavior for ui > 7 2 t> -1 is the behavior for 
usual 2-channel Kondo fixed point. 

Similarly we can also calculate the dot-spin correlation 
exactly. In the presence of the Ct term, the correlation for 
dot-spin has the usual 2-channel Kondo behavior 1/t for 
t < v/j 2 , and crosses over to 1/t 2 behavior for t > w/ 7 2 . 

We would like to repeat that all of the above results 
are valid only when 7 2 /k and T are much less than the 
Kondo temperature T2K ■ 

XGW is supported by NSF Grant No. DMR-97-14198 
and by NSF-MRSEC Grant No. DMR-94-00334. 



[2] 
[3] 

[4 

[5 
[6 

[7 

[8 
[9 

[10 

[11 
[12 

[13 
[14 

[is; 

[16 

[17; 

[is: 



D. Goldhaber-Gordon, et al, Nature, 391, 156 (1998); 

LANL archives cond-mat/9807233. 

S. M. Cronenwett, et al, Science, 281, 540, (1998). 

T. K. Ng and P. A. Lee, Phys. Rev. Lett., 61, 1768 

(1988). 

L. I. Glazman and M. E. Raikh, JETP Lett., 47, 452 
(1988). 

A. Kawabata, J. Phys. Soc. Jpn., 60, 3222 (1991). 

Y. Meir, et al, Phys. Rev. Lett., 66, 3048 (1991); 70, 

2601 (1993). 

H. Hershfield, et al, Phys. Rev. Lett., 67, 3720 (1991). 
A. L. Yeyati, et al, Phys. Rev. Lett., 71, 2991 (1993). 
N. S. Wingreen and Y. Meir, Phys. Rev. B, 49, 11040 
(1994). 

Y. Wan, et al, Phys. Rev. B, 51, 14782 (1995). 

I. Affleck JVucJ. Phys. B336, 517 (1990) . 

I. Affleck and A.W.W. Ludwig Nucl. Phys. B360, 641 
(1990) . 

J.M. Maldacena and A.W.W. Ludwig, Nucl. Phys. B506, 
565 (1997) 

V.J. Emery and S. Kivelson, Phys. Rev. B 46, 10812 
(1992) . 



Jinwu Ye, fcond-mat/9612029] ; Nucl. Phys. B512, 543 
(1998) 

P. Coleman, L. Ioffe and A.M. Tsvelik, Phys. Rev. B52 
6611 (1995). 

C. de C. Chamon, D.E. Freed, and X.-G. Wen, Phys. 

Rev. B 53, 4033 (1996) 

A. Lopatnikova and X.G. Wen, to appear. 



4 



